ELSEVIER 


Available  online  at  www.sciencedirect.com  - 

JOURNAL  OF 

V  ScienceDirect  POWER 

SOURCES 

Journal  of  Power  Sources  179  (2008)  424-429 

www.elsevier.com /locate /jpowsour 

Short  communication 

Simulation  of  a  thermal  battery  using  Phoenics® 


Giancarlo  C.S.  Freitas  a,  Fernando  C.  Peixotob,  Ardson  S.  Vianna  Jr. a’* 

a  Instituto  Militar  de  Engenharia,  Segao  de  Ensino  de  Engenharia  Qui'mica,  Praga  General  Tiburcio,  80, 

Rio  de  Janeiro,  Brazil* 1 

b  Universidade  Federal  Fluminense,  Escola  de  Engenharia,  Departamento  de  Engenharia  Qui'mica  e  de  Petroleo, 
Rua  Passo  da  Patria,  156  -  Bloco  D  -  Sala  307,  24210-240  Niteroi,  RJ,  Brazil 2 

Received  19  October  2007;  received  in  revised  form  12  November  2007;  accepted  13  November  2007 
Available  online  4  December  2007 


Abstract 

Thermal  batteries  are  primary  disposable  systems  specially  designed  to  develop  a  high  energy  density  in  a  short  period.  In  the  present  work, 
the  modeling  of  heat  generation  and  propagation  within  three  Ca/CaCr04  thermal  batteries  has  been  carried  out,  using  a  transient  model.  The 
commercial  CFD  software  Phoenics®  has  been  used  and,  through  a  typical  finite  volume  approach,  the  related  2D  transport  equations  have  been 
solved,  giving  the  time-dependent  temperature  profiles.  To  check  the  quality  of  the  model,  the  temperature  of  pseudo-equilibrium  state  (a  typical 
thermodynamic  parameter),  has  been  analyzed.  The  temperatures  related  in  the  literature  were  close  to  those  calculated.  The  results  also  indicate 
that  the  fusion  of  electrolytes  is  virtually  an  instantaneous  process  when  compared  to  the  time  to  reach  the  pseudo-equilibrium  state,  which  indicates 
that  the  generation  of  electrical  current  occurs  immediately  after  the  thermite  burning. 

©  2008  Published  by  Elsevier  B.V. 
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1.  Introduction 

Thermal  batteries  are  special  electrochemical  systems  in 
which  a  great  amount  of  energy  is  delivered  during  a  relatively 
short  period  of  time  [1].  The  cathode  and  the  anode  are  juxta¬ 
posed  to  an  electrolyte,  which,  at  room  temperature,  is  solid  and 
inert.  The  start  up  is  made  by  an  electrical  Squibb  that  initiates 
a  pyrotechnical  material  (thermite)  placed  between  the  cells.  As 
can  be  seen  in  Fig.  1 ,  the  internal  temperature  rises  quickly,  melt¬ 
ing  the  electrolyte  and  starting  electrochemical  reactions,  which 
generates  current  densities  much  greater  than  those  of  typical 
batteries  [2]. 

The  advances  in  this  technology  were  due  to  the  employment 
of  improved  quality  materials,  cathodes  and  anodes  in  particular 
[2] .  These  improvements  allowed  the  construction  of  more  com¬ 
plex  cells,  which  demanded  laboratory  tests  to  optimize  their 
performance.  The  costs  and  time  demands  in  assembling  and 
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testing  these  more  complex  batteries  are  considerable  high  and 
computational  simulations  prior  to  experimental  testing  are  of 
great  assistance  [3-6] . 

The  main  advantages  of  thermal  batteries  are:  long  shelf  life 
time  (up  to  25  years),  high  reliability  and  improved  resistance 
when  exposed  to  dynamic  charges,  such  as  shock,  vibration, 
acceleration  and/or  spin.  These  characteristics  allow  applica¬ 
tions  in  defense  devices  and  airspace  apparatus,  such  as  missile 
guidance  and  control  and  additional  energy  supply  for  airspace 
missions  and  satellites  [7]. 

The  thermal  battery  is  an  strategic  technology  mainly  used 
in  warfare  systems,  which  causes  a  low  availability  of  related 
scientific  information  in  the  open  literature,  most  of  the  infor¬ 
mation  being  obtained  from  limited  access  technical  reports, 
such  as  those  held  in  Sandia  Laboratory  [8,9],  QinetQ  (for¬ 
merly  Defence  Evaluation  and  Research  Agency — DERA)  [10] 
and  Brazilian  Army  Technological  Center  [2].  As  open  refer¬ 
ences,  we  can  mention  the  work  of  Nikitin  et  al.  [11],  which 
solves  the  heat  balance  equation  in  filter-press  type  batteries, 
Wu  et  al.  [12],  which  studied  the  heat  dissipation  of  lithium 
batteries,  Butler  et  al.  [13],  which  developed  refinements  use¬ 
ful  to  process  multiple-cell  and  multi-tap  testing  and  Schoeffert 
[14],  which  modeled  the  heat  transfer  within  a  Li/FeS2  thermal 
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Fig.  1.  Percentage  of  voltage  in  open  circuit  vs.  current  density  (mAcm  2). 


battery  and  compared  simulation  results  to  experimental  mea¬ 
surements. 

In  the  present  work,  a  transient  model  of  heat  generation 
and  propagation  within  a  Ca/CaCr04  (the  so-called  second  gen¬ 
eration)  thermal  battery  is  carried  out,  using  a  computational 
fluid  dynamics  (CFD)  code.  Briefly,  CFD  codes  are  computer 
programs  devoted  to  solve  the  systems  of  equations  governing 
heat,  mass  and  momentum  transfer.  These  systems  are,  typically, 
differential  balance  equations  coupled  with  empirical  rate  equa¬ 
tions  (for  example:  Newton’s,  Fick’s  and  Fourier’s  laws).  The 
software,  also,  must  use  a  suitable  numerical  method  for  the  task 
(for  example:  finite  element  or  finite  volume  methods)  and  pro¬ 
vide  an  adequate  graphical  interface  for  both  information  input 
and  output. 

The  simulations  indicate  that  maximum  electrolyte  tempera¬ 
tures  are  similar  to  the  temperatures  mentioned  in  the  literature  as 
those  produced  by  their  open  air  burning.  Besides,  transient  tem¬ 
perature  profiles  are  similar  to  the  results  published  by  Schoeffert 
[14]. 

2.  Fundamentals 

The  electrochemical  system  used  in  this  work  is  the  so-called 
second  generation  thermal  battery,  composed  of  a  Ca/CaCr04 
electrochemical  scheme.  It  works  at  an  average  temperature 
of  600  °C  and  its  main  elements  are  briefly  described  as  fol¬ 
lows. 

Thermite  is  the  heat  source  responsible  for  melting  the  ther¬ 
mal  battery  electrolyte,  leading  to  device  activation.  It  is  a 
mixture  of  zirconium  and  barium  chromate  sustained  by  asbestos 
fibers  in  a  circular  shape  [1,2]. 

The  X-ray  diffraction  analyses  of  ashes  resulting  from 
Zr/BaCr04  thermite  combustion  in  an  argon  atmosphere  sug¬ 
gests  a  hypothetic  decomposition  reaction  as  formulated  below 

[3]: 

3Zr  +  4BaCrC>4  3ZrC>2  +  4BaO  +  2&2O3 


With  this  equation,  thermal  effects  can  be  calculated,  since  the 
standard  heat  of  reaction  at  298  K  (A//|98)  is  [3] 

A//298  =  -476.3  kcal  mol-1  of  BaCr04 

The  cathode  is  made  by  applying  the  electrochemically  active 
material,  CaCr04  and  SiC>2  (also  known  as  cathodic  mixture) 
over  a  disc  of  pure  nickel.  To  achieve  this,  a  liquid  matrix  of 
trichloroethylene  and  dimethylphtalate  is  employed.  The  anode 
is  made  of  pure  metallic  calcium  and  the  electrolyte  of  the  ther¬ 
mal  battery  is  an  eutectic  mixture  of  inorganic  salts  applied  on 
a  glass  fiber  cloth  for  retention. 

3.  Methodology 

3.1.  Transient  model  construction 

Three  thermal  battery  designs  are  considered:  3-,  6-  and  12- 
element  cells.  In  a  3 -element  battery,  for  instance,  there  are  seven 
thermitic  elements  and  nine  electrolytes. 

The  transient  model  considered  involves  the  ignition  and 
the  heat  propagation  until  a  pseudo-equilibrium  temperature  is 
achieved. 

Each  battery  element  was  created  in  the  software  in  its  real 
dimension,  respecting  their  true  colors  and  using  the  materials  in 
the  software  library.  Thermites  were  defined  as  transient  sources 
of  energy  with  a  thermal  behavior  that  will  be  discussed  later. 
Additionally,  all  thermites  within  the  battery  were  assumed  to  be 
ignited  at  the  same  time  and  to  burn  at  the  same  rate.  The  fact  that 
all  elements  are  pressed  together  was  respected,  leaving  no  void 
spaces  between  them.  The  space  between  the  active  elements 
and  the  outer  case  was  “filled”  with  fiberfrax,  defined  as  the 
domain  material. 

The  natural  choice  of  the  coordinate  system  (. R ,  0  and  Z)  was 
respected  and  the  virtual  domain  had  the  dimension  of  the  steel 
case,  differing  only  in  height,  according  to  the  type  of  battery  in 
question  (3,  6  or  12  unitary  cells). 

The  heat  transfer  outside  the  battery  is  due  to  free  con¬ 
vection  and  radiation.  In  a  typical  battery,  in  which  diameter 
and  height  are  both  3  in.,  a  maximum  of  around  300  W  of 
heat  dissipation  by  (free)  convection  is  expected.  In  the  same 
way,  the  heat  dissipation  due  to  radiation  is,  coincidentally, 
around  300  W,  assuming,  in  both  cases,  that  the  surface  of 
the  apparatus  could  achieve  temperatures  as  high  as  500  °C. 
Since  a  single  thermite  generates  around  12,000  J  [3]  of  heat 
in  approximately  half  second,  the  system  can  be  considered 
adiabatic  for  our  purposes.  Additionally,  additional  thermites 
placed  at  the  top  and  the  bottom  exactly  compensate  this  heat 
loss  [14]. 

In  the  real  apparatus,  the  primary  ignition  is  due  to  squibs 
attached  to  wall  propagators  made  of  thermite  strips.  Therefore, 
ignition  of  all  the  thermite  happens  almost  at  the  same  time.  So, 
the  initial  simulation  time  was  considered  to  be  the  one  in  which 
all  the  thermites  are  ignited;  the  dynamics  of  the  propagators 
burning  and  their  ignition  by  the  squibs  are  outside  the  objectives 
of  this  work  (besides,  all  these  events  happen  in  a  very  small 
period  of  time). 
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High  temperature  variations  are  expected  in  the  thermal  bat¬ 
tery,  so  the  hypothesis  of  constant  thermal  conductivities  and 
heat  capacitances  does  not  apply  to  the  model  developed  here. 
As  a  result,  mean  values  for  those  entities  (for  each  substance) 
were  calculated  for  an  arbitrarily  chosen  temperature  interval 
from  298  K  to  902  K,  assuming  that  this  is  the  temperature  range 
in  the  adiabatic  system. 

The  burning  velocity  of  Zr/BaCr04  mixture  was  assumed  to 
be  0.56  ms-1;  attempts  to  accurately  measure  this  parameter 
were  tried,  but  until  the  end  of  this  work,  no  conclusive  results 
could  be  found.  As  there  are  two  propagators  in  the  battery 
and  they  act  on  two  diametrically  opposite  points  simultane¬ 
ously,  the  burning  time  can  be  calculated  as  t-  radius/burning 
rate  =  5.04  x  10-2  s. 

As  already  mentioned,  previous  works  [3]  employed  a  ther¬ 
modynamic  approach  to  calculate  the  heat  produced  by  a  single 
thermite,  which  was  found  to  be  12,140  J. 

The  heat  source  of  the  thermal  battery  was  modeled  as  a  pulse 
function,  which,  on  Phoenics®,  must  was  defined  as  shown  in 
Fig.  2. 

The  final  composition  of  each  substance  in  cathode  and 
anode  was  calculated  and  the  temperature-average  thermal  con¬ 
ductivity  and  heat  capacitance  were  obtained.  In  Phoenics® 
library,  two  new  materials  were  created  and  labeled  “cath¬ 
ode”  and  “anode”  with  their  respective  properties  calculated  as 
described. 

3.2.  Electrolyte  modeling 

The  value  for  the  KCl-LiCl  eutectic  mixture  heat  of  fusion 
was  the  one  mentioned  in  the  comprehensive  work  of  Masset 
and  Guidotti  [15].  Phase  transition  was  modeled  considering 
that  the  fusion  occurs  almost  instantaneously  (less  than  0.5  s) 
and  at  a  constant  temperature.  This  assumption  was  made, 
because: 

(a)  the  thinness  of  the  electrolyte  (0.10  mm)  allows  one  to  adopt 
the  hypothesis  that  this  element  fuses  completely — there  is 
no  need  to  establish  a  fusion  mushy  region  in  such  a  thin 
material; 

(b)  the  fusion  occurs  at  constant  temperature,  since  the  elec¬ 
trolyte  is  an  eutectic; 

(c)  the  heat  propagation  into  the  thermal  battery  occurs  mainly 
by  diffusion  and  convection  can  be  neglected. 


Additionally,  it  must  be  said  that  the  void  arising  from  the  sep¬ 
arator  processing  (electrolyte  +  binder)  for  the  Li/FeS2  reported 
in  [16]  is  minimum  in  the  present  battery,  and  corresponding 
thermal  effects  were  not  be  taken  into  account  in  the  present 
system. 

The  remaining  elements  of  the  thermal  battery  had  their  prop¬ 
erties  estimated  according  to  the  material  each  of  them  was  made 
of.  This  is  the  case  of  the  activation  base  (made  of  bakelite), 
mica  insulation  sheets,  separators  (made  of  steel),  negative  pole 
(a  nickel  disk)  and  also  fiberfrax  that  fills  all  the  empty  spaces 
after  device  mounting.  All  these  elements  could  be  found  in  the 
original  library  of  Phoenics®,  besides  fiberfrax,  the  properties 
of  which  were  entered  using  the  available  supplier  data-sheet. 

4.  Results  and  discussion 

The  results  of  the  3 -element  battery  simulation  will  be 
discussed  in  greater  detail  only  to  show  how  the  present  method¬ 
ology  can  be  applied.  Analogous  analysis  could  be  carried  out 
for  the  6-  and  12-element  batteries  but  we  opted,  for  sake  of 
simplicity,  to  merely  mention  the  results  for  both. 

The  simulation  of  the  3 -element  battery  showed  that  the 
pseudo-equilibrium  temperature  of  the  system  (PET)  is  around 
442  °C,  but  the  optimum  operation  temperature  (OOT)  of  this 
battery  was  found  to  be  between  550  °C  and  650  °C  according  to 
[2].  So,  to  ensure  enough  ionic  mobility,  the  3-element  thermal 
battery  could  have  its  number  of  thermites  increased  for  bet¬ 
ter  efficiency.  However,  according  to  a  more  recent  work  [15], 
which  mentioned  the  OOT  to  be  between  350  °C  and  550  °C, 
the  present  configuration  would  be  already  adequate. 

In  fact,  as  can  be  seen  in  these  simulations  -  and  could  even 
be  anticipated  by  the  device  anisotropy  and  by  the  existence  of 
refractory  materials  in  it  -  a  situation  of  thermal  equilibrium 
only  can  be  achieved  in  a  much  longer  time  than  that  for  the 
total  consumption  of  the  material  with  electrochemical  activity. 
The  OOT,  however,  was  defined  in  the  literature  [2,15],  so  as 
to  establish  a  parameter  for  calculation  of  the  best  amount  of 
thermites.  In  this  way,  its  comparison  with  the  values  of  pseudo¬ 
equilibrium  temperature  (PET)  obtained  in  this  work  will  serve 
a  qualitative  and  not  a  conclusive  purpose  for  the  acceptance  of 
the  experimental  arrangement  in  question. 

In  Figs.  3-7  it  can  be  observed  the  time  evolution  of  temper¬ 
ature  profiles  for  the  3 -element  battery,  respectively  at  instants 
4.44  s,  259.40  s,  745.50  s,  1 121.00  s  and  5505.00  s. 

The  simulation  revealed  that  temperatures  developed  next  to 
thermite  locations,  after  its  complete  burning,  is  much  greater 
than  the  equilibrium  temperature,  indicating  that  the  elements 
next  to  the  thermite  (for  instance,  electrolytes)  are  subject  to 
temperature  peaks  that  can  reach  high  values  (up  to  1000  °C). 
This  fact  is  due  to  the  high  combustion  heat  of  thermites,  asso¬ 
ciated  with  its  burning  rate  and  by  the  final  pressing  operation 
of  the  device.  Exhibiting  only  the  objects  of  interest,  Figs.  8-1 1 
better  illustrate  that  electrolytes  heat  quickly  after  the  thermite 
ignition,  allowing  the  time  required  for  the  eutectic  KCl-LiCl 
fusion  to  be  estimated. 

Temperature  profiles,  similar  to  the  ones  depicted  in 
Figs.  3-11,  were  found  for  the  6-  and  12-element  batteries  but 
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Fig.  3.  Temperature  profile  in  t= 4.44  s  for  the  3-element  battery.  Fig.  6.  Temperature  profile  in  t=  1121.00  s  for  the  3-element  battery. 
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Fig.  4.  Temperature  profile  in  t  =  259.40  s  for  the  3-element  battery. 


Fig.  7.  Temperature  profile  in  t  =  5505.00  s  for  the  3-element  battery. 


the  inclusion  of  the  corresponding  18  figures  were  considered 
unnecessary,  since  they  serve  almost  only  for  qualitative  pur¬ 
poses. 

It  can  be  noticed  that  in  t=  1.2  s,  the  electrolyte  temperature 
exceeds  the  fusion  temperature  of  eutectic  (350  °C).  This  shows 
that  that  the  activation  of  the  battery  is  almost  instantaneous, 
which  is  also  observed  in  our  laboratory  tests. 

The  total  computational  processing  time  was  around  6  h  in  an 
Atlhon  XP  1000  MHz  with  256  Mb  of  RAM  for  each  simulation 


presented  here.  Results  for  all  three  configurations  can  be  seen 
in  Table  1. 

From  Table  1,  it  can  be  seen  that  the  3 -element  battery  is 
the  only  one  which  presented  a  PET  outside  the  acceptance 
range  for  OOT.  This  fact  indicates  that  the  number  of  thermites 
for  such  a  device  should  be  increased,  in  order  to  improve  its 
electrochemical  behavior. 


Fig.  5.  Temperature  profile  in  t  =  745.50  s  for  the  3-element  battery. 


Fig.  8.  Electrolyte  temperature  profile  in  t  =  0.48  s  for  the  3-element  battery. 
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Fig.  10.  Electrolyte  temperature  profile  in  t  =  1.20  s  for  the  3-element  battery. 

To  evaluate  the  effect  of  increasing  the  numbers  of  ther¬ 
mites  for  this  device,  two  additional  simulations  were  conducted, 
adding,  respectively,  one  and  two  thermites,  similar  to  those 
already  present  in  the  battery.  The  analysis  times  were  the 


Fig.  11.  Electrolyte  temperature  profile  in  t  =  2.00  s  for  the  3-element  battery. 


Table  1 
OOT  x  PET 


Number  of  elements 

Calculated  PET  (°C) 

OOT  (°C) 

3 

440-442 

550-650 

6 

550-555 

550-650 

12 

643-684 

550-650 

Table  2 

OOT  vs.  PET  for  the  3-element  battery  with  supplementary  thermites 

Number  of  supplementary  thermites 

Calculated  PET  (°C) 

OOT  (°C) 

0 

440-442 

550-650 

1 

500-505 

550-650 

2 

550-556 

550-650 

same  previously  employed  and  the  results  are  presented  in 
Table  2. 

The  temperature  versus  time  profiles  for  electrolytes,  cath¬ 
odes  and  anodes  had  shown  no  appreciable  change,  so  the 
addition  of  one  thermite  may  raise  the  working  temperature  - 
and,  therefore,  improve  electrochemical  behavior  -  without  sub¬ 
mitting  electrolytes  or  electrodes  to  high  temperatures,  which 
could  cause  damage  or  undesired  parallel  reactions. 

On  adding  two  thermites,  however,  a  change  in  the  tempera¬ 
ture  profile  of  one  electrolyte  was  noticed,  which  heated  much 
faster.  Nevertheless,  the  achieved  PET  in  this  case  was  556  °C, 
within  the  acceptance  range,  which  indicates  an  improvement. 

Therefore,  simulations  show  that  it  is  possible  to  improve  3- 
element  battery  behavior  by  adding  one  or  two  thermites,  but  the 
advantages  and  disadvantages  of  each  case  must  be  taken  into 
account. 

When  adding  only  one  supplementary  thermite,  a  PET  below 
the  ideal  range  will  be  achieved,  but  the  risk  to  the  elements  next 
to  thermites  is  lower;  whereas  adding  two  thermites,  one  could 
raise  the  temperature  to  the  ideal  range  mentioned  in  the  litera¬ 
ture,  but  at  least  one  electrolyte  will  be  subject  to  an  undesired 
excessive  heating,  which  could  jeopardize  the  battery  perfor¬ 
mance. 

5.  Conclusions 

Despite  thermal  batteries  having  been  invented  in  the  1940s, 
their  improvement  and  optimization  is  still  the  subject  of  study 
and  research.  Since  it  is  a  device  specially  designed  for  defense 
applications,  scientific  articles  are  scarce,  incomplete  and  do  not 
reveal  technical  details  for  a  complete  technology  description. 

Concerning  model  construction,  the  results  in  this  work  indi¬ 
cate  that  the  simplifications  employed  and  the  methodology 
adopted  were  adequate,  because  simulation  results  replicated 
those  available  in  the  literature,  such  as  small  battery  action 
time — around  1.0  s  or  less.  Also,  the  maximum  electrolyte  tem¬ 
peratures  (those  next  to  thermites)  were  similar  to  thermite 
adiabatic  flame  temperature  (the  one  produced  by  open  air  burn¬ 
ing)  [3]. 

The  pseudo-equilibrium  temperature  of  Ca/CaCrC>4  batteries 
produced  by  the  Brazilian  Army  Technological  Center  (CTEx) 
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and  modeled  in  this  work  are  within  or  quite  close  to  the  opti¬ 
mum  operation  temperature  range  related  in  the  literature.  It  can 
be  concluded  that  the  dynamic  simulation  of  such  batteries  are 
similar  to  the  results  published  by  Schoeffert  [14]  taking  into 
consideration  the  particular  characteristic  of  each  system. 
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